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Abstract 

The vacuum stability lower bound on the mass of the Higgs boson is numerically 
investigated in an SU(2)l <8> SU(2)r symmetric Yukawa model, which describes two 
heavy degenerate fermion doublets in the limit of vanishing gauge couplings. Good 
agreement with perturbation theory is found, although the couplings are strong. 
The upper bound on the fermion mass and renormalized Yukawa coupling is also 



proven. 



O ! determined in the part of bare parameter space where reflection positivity has been 

m 

1 Introduction 

Cut-off dependent upper bounds on renormalized couplings arise in the Standard Model for non- 
asymptotically free couplings, if their continuum limit is trivial. Such "triviality bounds" show 
up in perturbation theory as an apparent inconsistency [|TJ, and have been intensively studied 
on the lattice by non-perturbative methods (for a recent review see ||). In Yukawa models 
there is also a lower limit on the Higgs mass, which is called "vacuum stability bound" || . On 
the lattice this can be understood as due to quantum effects inducing a positive renormalized 
quartic scalar coupling even if the bare quartic coupling takes its lowest possible value, namely 
zero [|]. 

In previous papers || the triviality upper bound on the Higgs mass has been investigated 
by non-perturbative methods in a Yukawa model with chiral SU(2)l <8> SU(2)r symmetry. The 
lattice formulation is based on the mirror fermion action |J with exact decoupling of the mirror 
fermions from the physical spectrum 0, (Another lattice formulation of the same continuum 
"target theory" is possible by using reduced staggered fermions 0.) Since the Hybrid Monte 
Carlo method JTOf is used, the minimum number of fermion doublets which can be numerically 



simulated, is two. In this context it is important that we are working in the limit of vanishing 
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gauge couplings, and hence by charge conjugation a mirror fermion doublet can be transformed 
to a fermion doublet. In fact, only the vanishing of the SU(3) co i our ® U(l)h yP ercharge gauge 
couplings is necessary. The chiral SU(2)l gauge coupling can be introduced, because SU(2) is 
pseudoreal flD], 0. 

In the present paper we continue the non-perturbative investigation of the same SU(2)l <S> 
SU(2)r symmetric Yukawa model as in refs. |T3|, In particular, we concentrate on the 
vacuum stability lower bound and on the upper limit of the renormalized Yukawa coupling in 
the region of positive scalar hopping parameter, where reflection positivity, implying unitarity 
in Minkowski space, can been proven Besides the physically relevant phase with broken 
symmetry, the renormalized Yukawa coupling is also computed in the symmetric phase, in order 
to check previous observations of very strong couplings there JTJ . 



2 Vacuum stability bound 

The lattice action and the definition of different renormalized physical quantities closely follow 
our previous papers on chiral Yukawa models with mirror fermions |14|, |], 13, |5|. The lattice 
action is a sum of the 0(4) (= SU(2)l ® SU(2)r) symmetric pure scalar part S v and fermionic 
part Si},: 



S Sm ~\~ S\f/ 



(1) 



tp x is the 2 (g) 2 matrix scalar field, and ty x = (if) x , x%) stands for the mirror pair of fermion 
doublet fields (usually ip is the fermion doublet and x the mirror fermion doublet). In the usual 
normalization conventions for numerical simulations we have 



S v = Yl \ o Tr (<Pt<Px) + A 



:Tr {<pZ<Px) - 1 



- « Tr (vi+^x) 



(XxlPx) + (ipxXx 



±4 _ _ 

K [i^X+UlM + {Xx+ftlvXx) + T ({Xx+^x) - {Xx^x) + {^x+ftXx) - {^xXx) 



=±1 



i^Rx^Lx) + {^Lx^x^Rx) + G x (XRxVxXLx) + (XLx<Px~XR 



ix) } 



(2) 



Here K is the fermion hopping parameter, r the Wilson-parameter, which will be fixed to r = 1 
in the numerical simulations, and the indices L, R denote, as usual, the chiral components of 
fermion fields. In this normalization the fermion-mirror-fermion mixing mass is /i^ x = 1 — 8rK. 

The fermionic part is given here for a single mirror pair of fermions. For the Hybrid 
Monte Carlo simulation the fermions have to be doubled by taking the adjoint of the fermion 
matrix for the new species. Taking the adjoint transforms fermions to mirror fermions and 
vice versa, but as noted before, without SU(3) co i OU r ® U(l)h yP ercharge gauge couplings they are 
equivalent to each other. 

The numerical simulations were performed on 6 3 -12 lattices at a bare scalar quartic coupling 
A = 10~ 6 . The small positive value of A was chosen in order to be sure about the convergence of 
the path integral. The Yukawa coupling G x was kept at zero, for exact decoupling of the mirror 
doublets ||. This allows to stay with the fermion hopping parameter K near its critical value 
at K = 1/8, as described in ref. ||. In the broken (FM) phase at the fixed values G$ = 0.25 
and Gv, = 0.30 the scalar hopping parameter was tuned to achieve a scalar mass of about 
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Table 1: The main renormalized quantities and the bare magnetization (er) = (\ip\) for several 
bare couplings and K-values. Points labelled by small letters are our data at A = 10~ 6 which 
were all obtained on 6 3 ■ 12. Points labelled by capital letters are the data for the upper bound 
at A = oo. For the points with double capital letters the lattice size was 8 3 • 16, whereas single 
capital letters denote data on 6 3 • 12. All data are collected from typically 10000 trajectories, 
except for the points labelled a, i, A, C where only about 5000 trajectories were run. As G x = 0, 
the renormalized coupling Gr x was always zero within errors and is not included here. 





G^, 


K 


<*> 


Vr 






9R 


Gr^ 




a 


0.25 


0.095 


0.251(7) 


0.196(11) 


0.79(7) 


0.271(8) 


42(11) 


1.38(4) 


1.42(20) 


b 


0.25 


0.099 


0.393(7) 


0.241(10) 


0.59(9) 


0.407(9) 


18(3) 


1.69(4) 


1.56(16) 


c 


0.25 


0.101 


0.525(6) 


0.304(9) 


0.57(3) 


0.59(3) 


10(2) 


1.68(10) 


1.53(40) 


d 


0.3 


0.090 


0.438(6) 


0.27(1) 


0.75(7) 


0.55(6) 


23(3) 


2.04(10) 




e 


0.3 


0.095 


0.754(4) 


0.417(11) 


0.76(7) 


0.91(6) 


10(1) 


2.19(11) 




f 


0.3 


0.100 


1.260(5) 


0.67(2) 


0.84(6) 


1.503(12) 


4.8(8) 


2.24(7) 




g 


0.62 


0.0 


0.424(3) 


0.29(2) 


1.23(9) 


1.23(3) 


53(11) 


4.2(1.2) 


4.5(4.0) 


h 


0.63 


0.0 


0.469(3) 


0.351(11) 


1.65(13) 


1.25(16) 


63(16) 


3.5(4) 




i 


0.65 


0.0 


0.523(3) 


0.34(2) 


1.39(14) 


1.6(3) 


53(18) 


4.6(4) 




A 


0.3 


0.30 


0.439(1) 


0.400(13) 


1.17(7) 


0.55(2) 


26(7) 


1.36(6) 




BB 


0.3 


0.27 


0.270(2) 


0.25(1) 


0.77(3) 


0.342(2) 


31(4) 


1.35(6) 


1.36(12) 


C 


0.6 


0.18 


0.3524(18) 


0.36(2) 


1.36(10) 


0.86(8) 


36(6) 


2.4(3) 




DD 


0.6 


0.18 


0.3389(13) 


0.339(16) 


1.31(7) 


0.86(11) 


38(4) 


2.5(3) 


2.3(2.4) 



m,R a ~ 0.6 — 0.8 and a not too small fermion mass, k was always kept to be non-negative, in 
order to be sure about reflection positivity, that is unitarity in Minkowski space The last 
points in the positivity region were fixed at K — 0, and then G$ was tuned to obtain reasonable 
masses. The results are summarized in table [1|, where also such points are included where the 
masses are not yet sufficiently tuned. 

Tuning the scalar and fermion mass is, of course, important in order to be close to the 
critical line separating the broken (FM) and symmetric (PM) phases [15], and at the same 
time avoid strong finite size effect. This prevents us from going on the 6 3 • 12 lattice to a very 
small Yukawa coupling G^ <C 1, because then the fermion mass becomes too small. At the 
strongest Yukawa coupling (at k = 0) the minimum of the scalar mass on our lattice is quite 
large (above 1), therefore we could not achieve the desired value rriR a ~ 0.6 — 0.8. This is 
similar to the behaviour observed at A = oo ||. More generally, in the investigated range of 
Gj/j the qualitative change of the renormalized quantities between 4 3 ■ 8 and 6 3 ■ 12 lattices is 
quite similar to the one at A = oo. This allows to choose the points with label c,d and g as 
optimal in the three sets of points for the lower bound in table |IJ These are shown in figs. 
1 and 2 together with the one-loop perturbative estimates of the vacuum stability lower limit 
at the given cut-offs. In these figures also the data at A = oo are included, which show the 
behaviour of the triviality upper bound for the renormalized quartic coupling. For completeness 
we display the results for A = oo in table [l] too, because more statistics has been collected for 
some points since publication of ||. 

The agreement between numerical simulation data and the one-loop perturbative estimates 



3 



is remarkably good. In particular, at the strongest couplings, within our errors, there is practi- 
cally no difference in the renormalized couplings between A ~ and A = oo. This means that 
the strong Yukawa coupling alone is able to induce the maximal possible renormalized quartic 
coupling at the given cut-off. 

Table [L] also includes results for the renormalized Yukawa coupling Gj^L, defined in terms of 
a three-point function (see B). At tree level, and moreover in the one- loop approximation, it 
coincides with Gr^. The numerical results show that both couplings are consistent with each 
other, even in those points where the errors are relatively large. 



3 Yukawa coupling 

One can see in table [I] that in the broken (FM) phase the renormalized Yukawa coupling Gr^ 
grows with the bare coupling G^ roughly linearly, up to quite strong values above the tree 
unitarity limit \[2tx ~ 2.5. In previous numerical simulations in the symmetric (PM) phase 



of the U(1)l ® U(1)r [14[] and SU(2)l ® SU(2)r |l3j symmetric Yukawa models very strong 
renormalized Yukawa couplings were observed, as well. We would like to see how strong the 
renormalized Yukawa couplings can be in the PM phase in the mirror- fermion decoupling limit. 
Since we do not want to go to the region with k < where reflection positivity could not be 
proven [|J, the maximal possible value of Gr^ should come from tuning along the k = 
axis in the PM phase. We tune G$ such that the physical scalar mass is around 0.7 on the 
6 3 ■ 12 lattice to avoid large finite size effects. We also run the same point on 8 3 • 12 and 8 3 ■ 16 
lattices to see how various quantities (especially the renormalized Yukawa couplings) change 
with different volumes. 

The renormalized Yukawa coupling matrix Gr = diag(G^, Gr x ) is defined as 

Gr = f R (k) Z~ 1/2 Z^ 1/2 f R (k) , (3) 

— 1/2 ~ 

where the three-point Green's function ($\l/\l/)o and are defined in [|TB|. TR^k) is the 

momentum-space renormalized fermionic two-point vertex function defined at the smallest mo- 
mentum k = (0, /c 4 ), fc 4 = tt/T. Near k = the behaviour of T_R(fc) is 

T R (k) ~ * 74 k A + Mr , k,= sinh = sin(^) , Mr = ( ° ^ j , (4) 

where fiR is the renormalized fermion mass. In previous work we did the leading-order ap- 
proximation by neglecting the term ^74^4 in the measurement since T is large. However, from 
data in |T3|] , we suspect that it might cause some small but visible effects on Gr. There- 



fore we decided to improve the measurement of Gr by including this leading-order correction. 
Also, we remove a previous inconsistency in normalization in the definition of Gr in the PM 
phase. Before, the convention was such that the full scalar propagator at zero momentum 
G(0) = J2x,y{4>Sx4>Sy) / L 3 T , after renormalization, was normalized to l/m R . However, this does 
not correspond to the natural normalization convention we used in the broken phase || . There 
the convention corresponds to the one in the PM phase such that G(0) is renormalized to 4/mJj. 
In order to have a fair comparison of the renormalized Yukawa couplings in both phases, we 
decide to switch to this new convention of normalization for Gr shown in (Q). This implies 
that data on Gr's in |Tj| should be scaled down by a factor of two in this new normalization 
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convention. (A similar inconsistency in the U(1)l®U(1)r symmetric model [|14| can be removed 
by dividing the values of Gr there by a factor of y/2.) 

One should notice that, as shown in ([3]), scalar quantities appearing in the definition of Gr 
should be the renormalized scalar mass m# and the wave function renormalization factor Zr 
defined at vanishing momentum However, the scalar mass we actually use in the mea- 

surement is the physical scalar mass obtained from a cosh fit of the scalar field correlation 
function along the time direction. The wave function renormalization factor Z3, which is de- 
fined through the residuum of the pole of the propagator, cannot be determined reliably from 
our statistics. Therefore we define the wave function renormalization factor Z^ in terms of the 
mass m<£ and the susceptibility (5(0) by means of 

Z^ = m\ (5(0) . (5) 

In general, uir and Zr is different from m$ and Z$. In a weakly interacting system (e.g. pure 
scalar A0 4 theory), they are very close to each other Ul6| . In our SU(2) mirror-fermion model, 
the bare Yukawa coupling G^ is large, therefore there is no guarantee that they are still close 
to each other. On the other hand, the correct measurement of tur and Zr is crucial to the 
measurement of Gr^ and Gr x , we therefore also measure itlr and Zr to see how they differ 
from rrirf, and Z§. Since we are on the finite lattice, we estimate Zr and rriR from 

Z R l = [G(k)- 1 - GiO)- 1 }/^ , m\ = Jg- , (6) 

where 

k=(0,h), h = ^, P = 4sin 2 (^). (7) 

The results of the numerical simulations in the symmetric phase are collected in table 0. 
From the table one can clearly see that Gr x is now zero within error bars as expected from 
the shift symmetry at G x = || [13]. This indicates that the term we used to neglect does, 



indeed, have some visible effect on the renormalized Yukawa couplings. Meanwhile, the natural 
definition of the renormalized couplings are taken at zero momentum since there is no infrared 
singularity in the PM phase. But Gr's we measured according to (|3j) are actually defined at 
k = (0, 0, 0, 7r/T), which will approach zero in the infinite T limit. If the inverse propagator had 
a non-negligible curvature near zero momentum, this would influence the determination of Zr 
and hence of other renormalized quantities. We therefore took measurements at two different 
values of T, i.e.: T = 12, 16, to see how we can get Gr's at zero momentum by extrapolation. 
Our data show that Gr^ and Gr x basically stay unchanged as T goes from 12 to 16. We 
therefore believe that curvature effects in the propagator are not large and that the data we 
have on Grs are, to a good approximation, the renormalized Yukawa couplings defined at zero 
momentum. 

The fluctuations of Zr and mjj are quite large, as expected. Within error bars, they agree 
with Z^ and m^. This shows that our measurement of Gr^ and Gr x using m</> and Z^ can be 
taken as a good approximation to the Grs defined in (|3]). 

As one can see from table [| on our lattices the maximal value of Gr^ in the region of non- 
negative scalar hopping parameter is more than twice the value of the tree unitarity bound. 
This supports our previous results in the symmetric phase on smaller lattices and less statistics 



Ti, [13j. With the present correct normalization the difference between the maximal value of 
Gr^p in the symmetric and broken phase is not dramatic, although the values are still larger in 
the symmetric phase (see tables [TJ-@). 
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Table 2: The main renormalized quantities at A = oo, G x = 0, G^ = 1.09, k = and 
K = 2/19 in the PM phase on various lattices. Data on 6 3 • 12 lattices are collected from 20000 
trajectories while on 8 3 • 12 and 8 3 • 16 we have about 6000 trajectories each. 



Size 




m R 


Zcf, 


Zr 


Hr 


Gr^ 


Gr x 


6 3 • 12 
8 3 • 12 
8 3 • 16 


0.69(3) 
0.65(7) 
0.55(3) 


0.60(17) 
0.56(11) 
0.51(11) 


1.52(5) 

1.84(15) 

1.71(33) 


1.11(38) 
1.36(57) 
1.44(73) 


0.664(5) 
0.654(9) 
0.63(1) 


5.84(16) 
6.04(37) 
5.87(54) 


-0.03(7) 
-0.12(17) 
0.07(9) 



4 1/N expansion 

The results of numerical studies of the phase diagram of the model at both A = oo and A = 10~ 6 
are reported in ||. The particular interest in the phase structure at small A arose after it was 
reported [[17] that models with naive fermions do not exhibit a ferrimagnetic (FI) phase and 
that instead a first order phase transition is observed. Clearly the absence of any such transition 
is of great importance for the study of bounds on the couplings within our model. 

Our Monte Carlo investigations show that at least the physically relevant phase transition 
from the PM to the FM phase is second order. In particular, the magnetization varies smoothly 
across the transition. 

The numerical analysis of the phase diagram was supplemented by a 1/Nf expansion in 
leading order, where Nf denotes the number of fermion-mirror-fermion doublet pairs. 

Our strategy was to calculate the one- loop effective potential to leading order in 1/Nf as a 
function of the fluctuation field 

a x = 4> 4x -s-(-l)^+^s, (8) 

where </>4 X is the fourth real component of the scalar field ip x , and s and s are the positions 
of the minimum of the effective potential at tree level with respect to ip x and the staggered 
scalar field <p x , respectively. After performing a Fourier transformation to momentum space, 
we finally obtain the effective potential V^fj[(x(0), ?(vr)] as a function of the fluctuation field at 
both the zero- and 47r-corners of the Brillouin zone. The qualitative features of the transitions 
from PM to FM (PM to AFM ]15[) can now be studied via the dependence of V e s on the field 
5(0) (the field a(n)). 

A first observation is that to leading order in 1/Nf the effective potential is a quadratic 
function of a(0) and 5(tt), and therefore we do not expect a first order phase transition which 
would rather require a quartic term, resulting in a double-well structure of the potential. 

For A = we obtain estimates for the critical value of the scalar hopping parameter k from 
the two gap equations 

- 16 + 2 - 87V, C?J ^ I - T _ fJ |L__ = o (for PM *-* FM), (9) 

16/t cr + 2 = (for PM^ AFM), (10) 

where = sin(g^), and the integral is taken over the Brillouin zone. It is therefore only the 
transition PM to FM which is affected by fermionic contributions in leading order. Its transition 
line bends down for increasing Yukawa couplings and finally intersects the straight transition 
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line from PM to AFM at k ct = —1/8. Obviously, in leading order of the 1/iV/ expansion the 
qualitative behaviour of the phase transition lines is similar to all our numerical studies [El [To| . 

According to the 1/Nf expansion the FI phase exists in leading order, since solutions to 
the minimum of Kff at non-zero values of both <r(0) and cr(vr) are found to exist. Furthermore, 
the position of the minimum of V e g varies smoothly across all phase transition lines, therefore 
suggesting the existence of second order phase transitions only. 

The findings from the 1/Nf expansion together with the Monte Carlo results are strong 
evidence that the allowed region of renormalized couplings can safely be studied at A ~ as 
well. 

5 Conclusions 

The main conclusion of the numerical simulations of heavy fermions in the SU(2)l <S> SU(2)r 
symmetric Yukawa model is that the estimates of the upper and lower bounds on renormalized 
couplings obtained in one-loop perturbation theory work well. All qualitative features of the 
one-loop /3-functions are supported, including the fixed point in the ratio of the fermion to 
scalar mass ^R^/m Ra (see fig. 2). The upper limit of the renormalized couplings (at the fixed 
ratio of gR/G 2 R ^) is provided by our requirement of reflection positivity. If this would not be 
imposed, the line towards the upper right corner of fig. 1 would be continued, but probably 
not much beyond our values. The reason is the appearance of the ferrimagnetic (FI) phase 
transition in the region of negative scalar hopping parameters, somewhat beyond the maximal 
bare Yukawa coupling we consider. The existence of the FI phase at small and large values 
of the bare quartic coupling A is implied by both numerical simulations and 1/N expansion. 
These conclusions are also supported in a recent numerical simulation of the same continuum 
"target theory" as ours, by using reduced staggered fermions [ fji| . 

Although the observed qualitative behaviour is certainly consistent with the one-loop per- 
turbative scenario implying the triviality of the continuum limit, one has to keep in mind that 
present simulations are done at relatively low cut-offs. In particular, the evolution of the cou- 
plings towards smaller values at decreasing cut-offs should be investigated in the future. At 
present the renormalized Yukawa coupling can have quite large values (see tables In the 

symmetric phase it reaches more than twice the tree level unitarity bound with both scalar and 
fermion masses about equal to 0.5 in lattice units. 
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Figure captions 



Fig. 1. The data for the upper and lower bounds on g R as a function of G 2 R ^ together with 
the perturbative estimates for scale ratios A/m Ra = 3 (solid curve) and A/m Ra = 4 (dotted 
curve). Open points denote the data for the lower bound (points c, d, g in table [I]), whereas full 
symbols are data for the upper bound. The 6 3 • 12-lattice is represented by triangles, whereas 
points on 8 3 • 16 are denoted by squares. 

Fig. 2. The mass ratio ^R^/m Ra as a function of G R ^ in comparison with one-loop 
perturbative estimates for m Ra = 0.75 (dotted curve), = 1 (full curve) and = 1.25 
(dashed curve). The solid horizontal line represents the fixed point at infinite cut-off. The 
explanation of symbols is as in Fig. 1. 
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Figure 1 



